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ABSTRACT 

We develop a theoretical framework to evaluate the energy spectrum, stationary states, and dielectric susceptibility of 
two Jaynes-Cummings systems coupled together by the overlap of their respective longitudinal field modes, and we solve 
and characterize the combined system for the case that the two atoms and two cavities share a single quantum of energy. 
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1. INTRODUCTION 

The Jaynes-Cummings (JC) systerrPEl comprises a two-level atom (2LA) coupled to a single optical resonator mode. 
Although Jaynes and Cummings emphasized the robustness of the semiclassical description as compared to the fully 
quantum model, manifestations of quantum field effects are now ubiquitous^' 7 leading to optimism that JC systems will 
soon be coupled together in one- or more-dimensional lattices* 5 thereby yielding novel condensed-matter phenomena.®^ 
£^ ' The first step to achieving coupled JC systems is to create and study a double- JC (DJC) system. 

We develop a theoretical framework for the DJC system by calculating its stationary states, energy spectrum and dielec- 
Ql trie susceptibility, and we show that this quadripartite system comprising two atoms and two field modes has fascinating 
features. The system is effectively characterized by two independent parameters: g for the coupling rate between the 
resonator and the single atom and k for the coherent photon hopping rate between the two resonators and proportional to 

■ the overlap of the two resonator field modes. By varying g and k, quite different features emerge from the DJC system, 
f — , Although we are interested in general properties of this system, our focus here is specifically on the case that the the DJC 

■ system shares precisely one quantum of energy. 

00 

The JC system has an atomic dipole with frequency lu r coupled to a cavity with frequency uj c via the atom's electric 

■ dipole between ground state |g) and |e). The dipole coupling frequency is g, and the atom-cavity decoupling frequency is 
A = u- d — uj c . In the rotating-wave approximation, the JC Hamiltonian is (ft = 1): 



2. DOUBLE JAYNES-CUMMINGS SYSTEM 



H sc = uj c I a f a + - ) + —u) a (r z + g (a^a- + aa + ) , (1) 



where and a are creation and annihilation operators for the cavity field and a± and a z with spec(<7r) = ±1 are spin 
operators for the atoms. The energy spectrum of the JC system is 



JC(O) = _A/o ,„JC( 



A/2, w± w = nuj c ± vW + A 2 /4 (2) 

for the ground state and the ?i th JC doublet, respectively. The stationary states are the ground state singlet |0g) and the 
excited state doublets |±) n given by 

|+)„ ± i |-)„ = e^- |„, g ) + e ±ie " |n - 1, e) , 6 n = ta n - 1 (2. 9 ^/A)/2. (3) 

Two neighboring cells (JC cavities) have overlapping evanescent mode functions resulting in an intercavity hopping 
rate n depicted schematically in Fig. [2] For a« and <r, the field annihilation operator and atomic electron energy lowering 
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Figure 1. Scheme for the system comprising two resonators each containing a two-level atom en and a bosonic mode en for i = 1, 2. 



operator for the i JC system (i = 1, 2), the DJC Hamiltonian is 



H = H[ 
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(4) 



with t/ the total number of quanta shared between the two atoms and the two cavity modes. Thus, the pure-state Hilbert 
space Jf? is a union of subspaces J4?("> with definite overall particle number v. 

Although we are interested in characterizing this system for all v, our focus in this paper is solely on studying the v = 1, 
which is intriguing in its own right. The v = 1 case provides an enticing simplification: the field modes can be treated as 
two-level systems (known as 'qubits' in quantum information parlance). As the system has atoms coupled only to the field 
modes and the field modes coupled to each other, for v = 1, the system corresponds to a chain of four qubits as shown in 
Fig.[2]with at most one qubit in the upper state. 

The spectrum and stationary states for the cases v = 0, 1 can be solved in closed form. Let = {|ni, Cx, C2)} 
be a basis for J^^ v > with rii denoting the number of photons in the i th mode and c^ = (e,g) the state of the i th atom. For 
the trivial case of no excitation in the system v = 0, we have 08^' = { | OgOg) }. For one excitation {y = 1) there are four 
basis states 

= {|0e0g) , |lg0g) , |0glg) , |0g0e)}, (5) 



In this basis 
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(6) 



with I the 4x4 identity matrix. The Hamiltonian matrix is not diagonal due to the presence of the coupling n. For 
independent cells, k = 0, and then the H^> matrix becomes a block diagonal with two 2x2 matrices each corresponding 
to a single cell. 

By diagonalizing the Hamiltonians for u = 0,1, the corresponding energy spectra are determined to be 



,(0) 



-A = 2^ JC (°\ 42 =oj c -±(A + e K )±^ 



(7) 



with e = ±, The spectrum o/°) is composed of a single nondegenerate state that is independent of g and k, whereas the 
spectrum a/ 1 ) is composed of doubly degenerate pairs for k = and is non-degenerate for k ^ 0. The coherent photon 
hopping effect (k ^ 0) does not lead to avoided crossings of the energies. Rather coherent hopping shifts the crossing 
point from A = to A = ±«. Only the atom-cavity coupling rate g results in avoided crossings. This avoided-crossing 
phenomenon for g and not for k is shown in Fig. [2] which shows the eigenvalues u/(± as a function of A for k 7^ 0, g = 
(left frame) and k ^ 0, g ^ (right frame). Evidently the energies cross at A = ±k for g = 0, and g 7^ leads to avoided 
crossing. 

For zero coupling k = 0, the first-level spectral term reduces to two degenerate eigenvalues 



w e± |re=0 
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Figure 2. The v = 1 energy spectrum vs detuning A illustrating avoided crossing effect to occur only for g ^ and at A = ±k: 

(a) 3 = 0, K = 2 and (b) g = 1, k = 2. 

with e irrelevant. Thus, the spectral values for k = can be understood as sums of spectral values for each of the two 
isolated JC systems. 

Now consider the v = 1 spectrum for the limited case of g = and let A = 0. In this case, a/ -* = 0, and w^t |g=o=A 
values are u) c (doubly-degenerate) and ui c ± n. These spectra correspond to the spectra for coupled harmonic oscillators as 
expected. Thus, strong inter-cavity mode coupling (large k) compared to atom-cavity coupling (g) is expected to make the 
DJC system behave nearly like coupled harmonic oscillators with a perturbation due to atom-cavity coupling. 

In the other limit we consider small k with strong coupling g. For v = 1 we have the spectral values 

«%* Ue ±g- T ± w (9) 

The shift ±g is due to vacuum Rabi splitting, and ere/2 is normal-mode splitting due to inter-cavity coupling. The next- 
order shift ±K 2 /8g is analogous to the ac Stark shift: the field of one cavity induces frequency pulling on the other cavity, 
which is detuned by g due to vacuum Rabi splitting. Thus, the strong-g, weak-K limit is equivalent to a weakly driven 
strong-coupling JC model. 

We now focus on the eigenstates of the system for v = 0,1. The v = case is trivial, composed of a singlet |0g0g) 
whose energy is independent of g and k. The eigenstates of the v = 1 Hamiltonian are 



|±,r 6 ) = W ±(|lg0g)-e|0glg)) + We ± (|0e0g)-e|0g0e)), (10) 



in which 



I = , W f= , (ID 

' 2 + 2 (r e T y/T+rl) W2 + 2 (r e T v 71 ^?) 



with 
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(12) 



2 5 

The states |±, r e ) in Eq. ( fTOb are like W states, i.e., a superposition of one excitation in each of the four degrees of freedom. 

In general, the states ( fTOb are non-maximally entangled states with unequal weighting of the field and atomic states in 
the superposition. The states are maximally entangled for uf — = 1/2. A close look at Eq. (fTTT i reveals that this could 
happen only for r £ = 0, which, according to (fT2l . only takes place for A = ±k. Thus, maximally entangled four-qubits 
states can be created at the thresholds A = — n and A — n. 



3. LINEAR SUSCEPTIBILITY 



The spectral properties and stationary states reveal the nature of the DJC system, but ultimately these features need to be 
observed experimentally. One way to observe these properties is to measure the dielectric susceptibility for a probe field 
directed through both cavities sequentially with measurement of the output field. The susceptibility is especially important 
for characterizing the DJC system for two reasons: susceptibility is experimentally meaningful on a macroscopic scale (i.e., 
without needing to manipulate individual elements of the system such as a single atom or cavity), and the susceptibility 
would provide a direct signature of a quantum phase transitions in the JC lattice case. 

The connection between susceptibility and spectrum is as follows. Given that the system is prepared in some sector of v 
excitations, the probe field will excite the system only if the field frequency is close to resonant with a transition from this 
energy state to a higher energy state in the system. The susceptibility is proportional to the probability of the transitions 
from v = 1 states to the ground state. For the system prepared in the sector v = 1, the susceptibility is given by 



xo 1 ' (wp) = E E r ^ ( u * -^p- ^) > w 



»=± e=± 

where 7 a is a small parameter describing a finite width of the transitions and 



r iie = 7 1 {i, r e \ (<r+ + a+) |0) | 2 + 7 c |<», r e \ (a\ + a|) |0) 



(14) 



is the total probability of the transitions from the energy states |±, r e ) to the ground state |0) = |0g0g) . The total probability 
of the transitions T^ is a sum of the squares of the absolute values of the amplitudes of transitions from the energy states 
|±, r e ) to the ground state |0) caused by spontaneous emission from the atoms, occurring with the rate 7, and by damping 
of the modes of the JC cavities with the rate 7c . 

The atomic dipole operators and the field operators have matrix elements 

(±,r e | {or\ + &l) |0) = (1 - e)wf, (±, r«| (a\ + a|) |0) = (1 - e)uf, (15) 
which cause transitions from the v = 1 states |±, r c ) to the ground state |0) to occur with probabilities 

r ± , € = (l-e) 2 (j\wt\ 2 +lc\ut\ 2 ). (16) 

Clearly the transition rates from states with e = +1 are zero so that the states are dark (non-radiative) states irrespective of 
the spontaneous emission, cavity damping, and whether the states are maximally entangled or not. This property contrasts 
with the two-qubit case for which a dark state can be created between the qubits only if the qubits are identical.^ 

In our case of four qubits, the atoms are degenerate in frequency and the field modes are also degenerate in frequency, 
but these frequencies differ from each other. Of course, the states e = +1 would not be dark states when we unbalance the 
symmetry between the atoms and/or the cavity modes, i.e. when either the atoms or the field modes are not degenerate in 
frequency. The symmetry could also be broken by allowing the atoms and/or the cavity modes to be damped with different 
rates. If the atoms are damped with different rates, say 71 and 72, and the cavity modes are damped with rates j c i and 7 C 2, 
respectively, then the transitions between the states |±, r c ) and |0) occur with probabilities 

r±, e = (VtT - \wf | 2 + (Vt^I - ^Vt^) 2 \ u t I • ( 17 ) 

Clearly, the transition probabilities are different from zero irrespective of e. In this case, the absorption spectrum is 
composed of four peaks of different magnitudes. 

For the case that 71 = j c i and 72 = 7 C 2, the probabilities ( TPTI i reduce to 

r±,e = (V^h-e^f. (18) 

In this case, transitions from the two states corresponding to e = +1 occur with the same probability (-1/71 — ^/72) 2 ■ 
Similarly, transitions from the states corresponding to e = — 1 occur with probability (^/tT + Jyi) 1 ■ Thus, the resulting 



absorption spectrum of a probe field is expected to be symmetric regardless of whether the states |±, r e ) are maximally 
entangled or not. Thus, if the atom and the cavity mode of cell 1 are damped with the same rates and also the atom and 
the cavity mode of the cell 2 are damped with the same rates, which may or may not equal the damping rates of cell 1, the 
transition probabilities are independent of whether the states |±, r e ) are maximally entangled or not. 

However, in the case of 7 7^ 7 C , the spectrum could be symmetric only if the states |±,r £ ) are maximally entangled. 
Otherwise, the spectrum is asymmetric. From Eqs. <TT~3T > and (116b . the imaginary part of the susceptibility (absorption 
spectrum) is a sum of two Lorentzians: 



1m 



7 a r_,- 



,(1) 



'7a 



(19) 



+ 7a 2 



The symmetry of the spectrum depends on the ratio between T + and r_ , and the spectrum could be symmetric only 
for T + ,_ = r_ t _. From Eq. dT6b we see that = T_ ; _ holds only if \uf | 2 = \wf f = 1/4. According to Eq. (fTBl , 
it happens only when the states |±, r_) are maximally entangled. Hence, an observation of the symmetric absorption 
spectrum when 7 7^ 7 C could be regarded as an indication of the presence of maximally entangled states in the system. 
Figure[3]shows the absorption spectrum for 7 7^ j c , k = 2g, and different A. We see that as long as A 7^ k, the spectrum is 
composed of two peaks of unequal amplitudes. The spectrum becomes symmetric at A = k. In this case the states |±, r_) 
are maximally entangled states. Thus, the symmetry of the spectrum can be regarded as an indication of the presence of 
maximally entangled states. Figure|4]shows the absorption spectra for 71 7^72and7 c i 7^ 7 C 2 with (71, 72) 7^ (7ci,7c2)- In 
this case the spectrum is composed of four peaks and is always asymmetric. However, one can notice from the figure that, 
at A = —k, two peaks corresponding to transitions from the states |±, r+) have equal amplitudes and are symmetrically 
located about (cu c — oj p ) / g = 0. Similarly, at A = K, two other peaks that correspond to transitions from the states |±, r + ) 
now have equal amplitudes and are symmetrically located about (cj c — uj p ) / g = 0. Again, equal amplitudes of the peaks 
indicate that the states which the transitions correspond to are maximally entangled states. 

The presence of the threshold values for A at which maximally entangled four-qubit states are created could be predicted 
from the structure of the Hamiltonian of the system. Instead, working in terms of the two coupled JC systems, we can 
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Figure 3. Absorption spectra Im j^Xo 1 ' (^p)j °f two coupled JC cells plotted as a function of (uj c — co p )/g for j a /g = 0.05, j/g = 0.01, 
7c/g = 0.02, n/g = 2 and different A: (a) A/g = 0, (b) A/g = -5, (c) A/g = 2 and (d) A/g = 5. 
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Figure 4. Absorption spectra Im ^Xo^ ( a 
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of two coupled JC cells plotted as a function of (u c — w p )/p for 7a/<? = 0.05, 
(0.2, 0.01), n/g = 2 and different A: (a) A/g = 0, (b) A/g = -2, (c) A/g = 2 and 



introduce symmetric and antisymmetric combinations of the atomic and field operators (collective modes) 

A\ = —7= {ax + a 2 ) , A 2 = -^= (&x - a 2 ) , 
V2 v2 

& = (K + ^2) > ^2 = ^= (of - 5-2 ) , (20) 
and find that, for any v, the Hamiltonian of the system can be written as 

H = {uj q + k)A{Ax + (uj c - k)A\A 2 + -U) & S Z + g (A\Sx + A\S 2 + H.c.) . (21) 

Thus, the DJC system is equivalent to two independent and non-degenerate collective systems corresponding to symmetric 
and antisymmetric combinations of the modes. The collective bosonic modes are coupled to the collective atomic system 
with the same coupling strength g. 

Hamiltonian ( |2TT i can be written as 

H = uj a (A\Ax + A\A 2 + ^sA - (A - k)A{Ai - (A + n)A\A 2 + g (A\Sx + A\S 2 + H.c.) . (22) 

with the first part representing the average free energy of the field modes and the atoms, the second part representing the 
shift of the energies of the superposition modes from the average energy, and the last part representing the interaction of 
the field modes with the collective atomic systems. 

Notice that, at A = re, the systems Ax and Si are resonant, so that the coupled two-qubit system can then be maximally 
entangled. Similarly, at A = —re, the systems A 2 and S 2 are resonant and therefore can be maximally entangled at that 
frequency. Clearly Eq. ( fTOb implies that, at A = —re, the avoided crossing occurs between states |±,r+) that involve 
antisymmetric combinations of the atomic and field states, whereas, at A = re, the avoided crossing occurs between 
states |±,r_). Thus maximally entangled states can be created in the system for A = re and A = — re. 



4. CONCLUSIONS 



We have constructed a framework for calculating the energy spectrum, stationary states, and dielectric susceptibility of 
two Jaynes-Cummings systems coupled together by the overlap of their respective longitudinal field modes and solved it 
for v = and v = 1 excitations of the system, which can be understood in terms of four coupled qubits. For weak coupling, 
the pair of systems is similar to a single Jaynes-Cummings system undergoing an AC Stark effect, and for strong coupling 
the behavior is similar to two coupled harmonic oscillators. For moderate coupling strengths, the pair of atoms and the 
pair of field modes can be highly entangled states, and, where the spectrum exhibits avoided crossings as a function of the 
detuning, the atoms and fields are found in maximally entangled four-qubit W-like states. We also show the susceptibility 
and the absorption of the system that explore the entangled features of the system. 
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